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Oh: 

^ ' Abstract 

(N 

^ ' Based on the infrared improved soft-wall AdS/QCD model which can lead to a consistent prediction 

OO \ for the mass spectra of light resonance mesons, we extend it to the finite temperature system and carry 

CO ■ out the calculation for the quark number susceptibility with finite quark mass. As a consequence, 

t^^ ' we show that the quark number susceptibility grows rapidly with a continuous blow-up in a narrow 

Ch ^ temperature interval as the temperature increases, and there is a peak when the temperature is around 

r~^ ■ T ~ 160 ~ 190 MeV depending slightly on the models, the resulting critical temperature is about 

f^ ' 170MeV, which agrees remarkably with the lattice QCD calculation. 



X 



1 Introduction 

Quantum chromodynamics describes the dynamics of quarks. Its gauge group is non-Abelian and the 
gauged gluon fields have self-interactions. As a result, the beta-function at ultraviolet (UV) region becomes 
negative that leads to the asymptotic freedom [1,2]. While at low energies, the non-perturbative dynamics 
such as chiral symmetry breaking and linear confinement appear as the gauge coupling constant grows. 

The duality between gravities and gauge theories [3-5] gives a new approach to the study of the low 
energy QCD. Among the studies of the duality, the ones developed phenomenologically are referred as 
bottom-up approach (AdS/QCD model) [6,7] (There are also good reviews, e.g. [8]). One of the important 
achievements in the AdS/QCD model are the realization of the chiral symmetry breaking in the hard- 
wall AdS/QCD model [6] and the linear confinement in the soft-wall AdS/QCD model [9]. But in the 
hard-wall model, the resulting excited meson mass spectra are found to be contrary to the experiments. 
In the soft-wall model, one can obtain desired mass spectra for the excited vector mesons, but the chiral 
symmetry breaking occurs inconsistently. Several groups studied on these problems [10,11]. A predictive 
soft-wall AdS/QCD model was constructed by simply modifying the 5D AdS metric in the infrared (IR) 
region [13] and the resulting mass spectra for the light resonance states of pseudoscalar, scalar, vector and 
axial- vector mesons remarkably agree with the experimental data. It is then interesting to study the phase 
transition for such a holographic AdS/QCD model with finite temperature, which comes to the purpose of 
our present paper. Thus we will extend such an infrared improved soft-wall AdS/QCD model [13] to a finite 
temperature system by considering a black hole metric and investigate the phase transition by computing 
the quark number susceptibility (QNS). The QNS is known to be the response of quark number density 
to infinitesimal change of the quark chemical potential as Xq = duq/dfiq. Using the fluctuation-dissipation 
theorem, it can be expressed as [14] 

X, = -limReGf,(c. = 0,fc), (1.1) 

where Gf((aj,fc) is the two-point retarded Green function for time-component of quark number current. 
From the lattice QCD studies [15,16], it strongly suggests that the behavior of Xq/T"^ on the transition 
between the hadronic phase and the Quark-Gluon-Plasma phase at zero quark chemical potential is a 
continuous but rapid in a narrow temperature interval around T ~ 170 MeV [17]. In this paper, we are 
going to compute the QNS with the constituent quark mass friq due to the dynamically generated chiral 
condensation of quark a ~< qq >, which is defined as fhq = rn^ + <t//U^ with niq the current quark mass 
and cr/ ^■^ regarded as the dynamically generated quark mass for the typical QCD scale /i r^ ^qcd due to 
QCD confinement [18], and to show how the behavior of Xq/T"^ changes as the temperature and dynamical 
quark mass and what is the critical temperature at zero chemical potential. 



This paper is organized as follows: in Sectl2l we present an introduction to our working model and 
the C/(l) gauge field dual to the quark number current. In Sect (31 we show how to get the quark number 
susceptibility from the solution of the C/(l) gauge field, and provide the numerical results with a detailed 
discussion. Our summary and remarks are presented in Sect|4l 

2 The models 

As the predictive soft-wall AdS/QCD model constructed by simply modifying the 5D AdS metric in the 
infrared (IR) region [13] can result the consistent mass spectra for the light resonance states of pseudoscalar, 
scalar, vector and axial-vector mesons, we shall extent such a model to a finite temperature system for 
studying the phase transition. This may simply be realized by considering a 5D black hole metric as follows 

ds' - a\z) (f{z)dt' + Y^dx^-j^Y (2.1) 

with 

a^{z) = l/z^ + ^l and f{z) = 1 - {z/zq)*, (2.2) 

where z has the relation with the normal radial coordinate r as z = l/r and zq is the horizon. The scale 
factor ^g characterizes the confinement of QCD in the infrared region by breaking the conformal symmetry, 
it is fixed by fitting with a linear function for the excited states and has been given in [13] from a consistent 
prediction for the resonance messon mass spectra. Hawking temperature is given as T = l/(zo7r). 

The action with the background dilaton ^{z) which has a nontrivial solution for forming the soft- wall 
is given as [13] 

S^^ J d'x^e-^^'^Tr{\DMX{x,z)\'-m^x\Xix,z)\^y (2.3) 

where m^ = —3 so that the solution of dilaton $(2;) obtained below does not have the singular at UV-limit 
(z — > 0). X{x, z) is the bulk scalar field, and we denote its vacuum expectation value (VEV) as X{z) for 
simplicity 

X{z)^^v{z)l2, (2.4) 

where I2 denotes 2x2 unit matrix. For our purpose in the present note, we consider that the VEV X{z) 
is taken as the same form given in [13], which is listed in Tabled) While, the dilaton $(2:) is fixed from 
the equation of motion for X(z) which is given as 







d, (a3(z)/(z)e-*(")a,z;(z)) - a^{z)e-'^^'^mj^viz). (2.5) 



where the resulting $(2:) from eq. (|2.5p will be different from the one yielded in [13] due to the inclusion of 
the black hole metric f{z). 

Model v{z) Parameters 

la z{A + Bz^){l + Cz^)-^ B^a/C + mqQC, C =^ B/{^i,a) 

lb z{A + Bz^)il + Cz^)-^/^ B^a/C+^mgCC C ^ {B^ / {^idl^)f''' 

Ila z(A + Bz2)(l + Cz4)-i/2 B = a/C, C = {B / {iidi)f 

lib z{A + Bz^){l + Cz^)-^/^ B = a/C, C = [B^ /{^ldl^)f''^ 

Table 1: A = iriq^ and ( = ^3/(2%) for all the models, while /i^ = 3/(2a)/ig with a = 1 for the la, Ila and a = 1/2 
for the lb, lib. A, B and C satisfy ea. (|2.7|l . Numerical values of mq,"/ and a are given in Table[2] 

Such an VEV may be regarded as the one at the low temperature phase. It was shown that the boundary 
behaviors of the VEV play an essential role [13] for mass spectra of resonances mesons. In general, its UV 
boundary behavior is fixed by the AdS/CFT correspondence, while its IR boundary behavior is imposed 
so as to obtain the required boundary behavior of dilaton at the IR- region. From Tabid, we can sec that 
v(z) behaves in the IR-limit z — > 00 and the UV-limit z — > as 

^WL^oo=T Y^Ms^ =7(Md^r, (2.6) 



a 



^Wluv-ihnit=™^C^ + 7^'- (2-7) 



where /i^ is defined with a = I for the Model la and Ila and a = 1/2 for the Model lb and lib, and ( is fixed 
as C = v3/(27r) [13]. niq and a denote the current quark mass and the chiral condensation respectively, 
and the constituent quark mass in this paper is defined as fhq = mq + a j ^i} [18]. Their numerical values 
are given in Tablel2l which have been fixed by minimizing the breaking of the Gell-Mann-Oakes-Renner 
relation as 

!l ml = 2mq a (2.8) 

at the 1% level and the experimental values: to^ = 139.6 MeV and f.^ = 92.4 McV. 

To compute QNS by using cq. p.ip . let us now introduce an U{\) gauge field which is dual to the 4D 
quark number current. Changing the coordinate to u = [z/ zh)"^ , and considering the following action 



S^-^J dud^X^e-^^'^^FranF'^^. 



(2.9) 



we can simply get the equation of motion = df^{y/ge "^F^") for A^. Performing the Fourier transfor- 
mation as Ap_{t,x,u) = /755TT e~""*''^^^''Afj_{uj,k,u), choosing the momenta k along the x-direction and 



Model la lb Ila lib 

m, (MeV) 4.16 4.64 4.44 4.07 

cr3 (MeV) 275 265 265 272 

7 0.178 0.136 0.153 0.112 

Table 2: ^ij = 3/{2a)^l = 445^ with a = 1 for the la, Ila and a = 1/2 for the lb, lib. 

taking the axial gauge Au{u) = 0, we obtain the following equations of motion for the needed components 

At : 0= (e"*y^5™'5"^0' - g'e-*V5.9"V*^t " coqe-^y^g^^g'^A,, (2.10) 

A„ : = wg"A;-95--A;, (2.11) 

where q and uj are defined as q = k/ (2'kT) and uj = u)/(27rT), the prime denotes derivative with regard to 
u, and in this study we use x-direction as a usual z-direction in four dimensional space-time due to what 
we have already used the symbol z as the fifth dimension here. In the above equations, we represent A^ 
as A^, and will apply this manner in this paper for simplicity. 

3 The results 

It turns out that the solution $(u) can be obtained analytically from the equation of motion eq. ()2.5p for 
all the modelcl. In obtaining the solution for the dilaton $(«), we take the integral constant so that the 
solution at the boundary behaves as 

*(^)luv-mnit = 3M^^ (3.1) 

Here we abbreviate to represent its explicit results due to a long expression. By using the solutions $(u) 
in all the models and the equation of motion cq. (|2.10|) . we are able to get the retarded green function. 
From eq. (|2.10p . we can obtain the following equation 

with ij} = e~*Y^(/""(7**^J. Note that there is a singular at u = 1 due to the horizon, and the deepest 
one is with the power —2. To get rid of such a singular, we assume that the solution of ^ has the form 
^ = (1 — uY'^'^ X{u^Lo^q) with the regular function X{u,Lj,q). Then by appropriately determining the 
index v, we can regularize the singular behavior. In general, there are two options for the values of t^, 
which corresponds to the in-going and out-going solutions respectively. In our present consideration, as 



*We use mathcniatica i 



the value of ui is taken to be zero eventually following eg. ljl.ip . it is not important to specify which one we 
are actually taking. In the following treatments, we consider the case that a;,*? ^ 1 in the hydrodynamic 
expansion. With such an assumption, the form of the function X{u^Lo^q) can be expressed in the general 
form: X{u, u, q) = i/jq + q^ijji{u) + uj'^ip2{u) + 0{uj'^, q'^). As a result, provided that i' ^- for w — > 0, we 
can show that -00 and "ipi satisfy the following equations 

^gXXgt 

v4 y t^ 



0=f ^ V '^ (3 4) 



It is then not difficult to find out the solution for -00 

V'o = Ci + C2 / due-'^^g^^g'K (3.5) 



Ji 
as the integrand ^/gg^^g^^ has the singularity arising from the horizon, we then take C2 = 0. And the 
solution for ipi is given 

^i = C,e-^^gg-g^' (^" d,'__l__ + z?i) , (3.6) 

From the condition tl''{u = 1) = 0, the integral constant Di can be fixed. From the Aj— component of 
ea. (|2.10p and the definition of ^jj with w — > 0, we have 



•^-^■(r^«'Fwv^+''')+°'"^''''^ 



(3J) 



where Ci can be fixed by using Km At = A\ to be 



C. ^ 4") [[ '-' e-^Jgg-9^^ + i..) ' + 0(.^ <Z^). (3.8) 

Here, Ai plays the role as the source for the 4D quark number current mentioned in the above eq. (|2.9p . 
From the definition of ■(/;, we can obtain the following expression 

^i = e-^^gguugu + ^(-'' '^') + ^(-'' '?')■ (3-9) 

We now come to discuss the retarded green function through the GKP-W relation [4, 5] and the 
prescription [12]. It turns out that the on-shcU action takes 



Thus the retarded green function Gf ^ is found to be 



So ~ ^, J (^ [^AA', - jA^A'^ 



(3.10) 



G^(m„^,,a) = -^ (/°^ \-^ Jg.ugU +^1) +Oi^',q% (3.11) 

5 



where a), is normalized to 1. By taking g — > with oj = 0, we can obtain the quark number susceptibihty 
Xq through eq. (jl.l[) . In the following computation, we will neglect the high order correction part of 
0(0;^, q^), and also take 95 to be gs = 1 for simplicity as the physics property of the green function Gf^ 
is irrelevant to its absolute value. 

We shall compute Xq/T'^ as function of three parameters a, nig and /i^. Typically, we will discuss the 
variation of a with a range ±25 toward the value of <j^'^ in Tablel^l This does not mean the uncertainty 
but the range to see how the QNS Xq/T"^ changes as the dynamical quark mass. On the other hand, we 
shall consider the uncertainty of about 10% toward the values of /i^ in Tablej^J which is mainly caused 
from the difference between the predicted mass spectra and experimental data for the resonance meson 
states [13]. 



la : (j^/^ = 275 ± 25 (m, ~ 4.16 T 1-2), 7 = 0.178 and ^id ^ 445 {fig 

lb : fji/^ == 265 ± 25 (to, ~ 4.64 =F 1-2), 7 == 0.178 and fid = 445 (fig 

Ila : cri/3 = 265 ± 25 (m, ~ 4.44 T 1-2), 7 = 0.178 and ^a = 445 {^ig 

lib : a^/^ = 272 ± 25 (m, ~ 4.07 T 1-2), 7 = 0.178 and fid = 445 {fig 



363) 

257) 
363) 

257) 



Table 3: The parameters used in Fig[21 7 is dimensionless and the unit of the rest are MeV. rUq is automatically 
fixed by the relation. (|2.8p . 



la 

lb 

Ila 

lib 



iq = 4.16, cri/3 = 275, 7 = 0.178 and fid = 445 ± 50 (fig ~ 363 T 40) 
iq = 4.64, cri/3 ^ 265, 7 = 0.136 and fid = 445 ± 50 {fig ~ 257 T 30) 
2g = 4.44, cri/3 ^ 265, 7 == 0.153 and fid = 445 ± 50 (fig ~ 363 T 40) 
iq = 4.07, cri/3 == 272, 7 = 0.112 and fid = 445 ± 50 {fig ~ 257 T 30) 



Table 4: The parameters used in Fig[3] 7 is dimensionless and the unit of the rest are MeV. The relation for fi^ 
and ^g are given in eq. 



Before making the computation of quark number susceptibility, we may check that the solutions A[ are 
well-defined in the entire bulk (FiglT|). 




Figure 1: The figures show the regularity of the solution A't in the entire bulk. Note the scale taken in each 
figures. The x-axis and y-axis denote the radial coordinate u and the real part of the solution A[ in the la model, 
respectively. On the other hand, the imaginary-part in all models always takes the value extremely close to zero. 
Considering it as just a numerical error, we then skip to show it here. The three figures are for the real-part of 
the model la: the left one indicates no-UV divergence, the center one shows the behavior in the whole region, and 
the right one shows finiteness and the boundary condition is taken correctly around the horizon. The behaviors 
for other models are the same. In the figures, we have taken u £ [10~ "", 1 ~ 10^^'] with the in-going boundary 
condition, and q and uj are taken to be zero. The parameters used here are given in Table[4] with temperature 
changing as T=120 (green, dashed), T=200 (red, dot-dashed), T=300 (blue, dotted), T=400 (magenta, solid), 
(Unit is MeV). 



From the results shown in Figll] it is seen that there is no divergence involved in our present considerations, 
therefore we do not need to care the regularization for the retarded green function [19,20]. Thus by using 
the parameters in Tablel3]and|4l we obtain the results shown in Figl2]and[3l 
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Figure 2: Results with fixing Hd and varying ruq and cr^ " as in TableO The different curves from top to bottom 
in all figures correspond to a''^^^ = 275 -I- 25 (green, dashed), cr^''^ = 275 (red, solid) and a''^^^ = 275 - 25 (blue, 
dot-dashed). The changes of a is just to see the dependence of QNS on the constituent quark mass given by 
■fhq = niq + a/fi^. It is observed that the appearing of peaks and its locations in all models agree with the results 
of lattice QCD [16]. 



It can be seen from Figsl2] and [3] that the QNS Xq/T'^ emerges rapidly with a continuous blow-up 
in a narrow temperature interval as the temperature increases. There is a peak when the temperature 
reaches to be around T ^ 160 ^ 190 MeV depending slightly on models, and after that, the behaviors of 
Xq/T^ get slowly in all the models. As mentioned in the introduction, the results of lattice QCD [15, 16] 
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Figure 3: Results with fixing niq and a^' as in Table[4]and varying ^d- The difTerent curves from top to bottom 
in all figures correspond to ^d = 445 — 50(green, dashed), /i^ = 445(red, solid) and /i^ = 445 + 50(blue, dot-dashed). 
The possible differences with regard to ^d arise from the one between the predicted resonance meson mass spectra 
and the experimental data. It can be seen that the appearing of peaks and its locations in all models are in good 
agreement with the lattice QCD' results. [16]. 



have shown that when the phase changes from the hadron phase to the quark-gluon-plasma phase at zero 
chemical potential, Xg/T^ grows rapidly and continuously with a peak around T ~ 170 McV. Therefore, 
the appearing of peaks and its locations in all models are in good agreement with the lattice QCD' results 
Let us now make a determination for the critical temperature in each model, which is simply fixed at 
the location of the peak 



la: 



173.6 ±3.2 1,, ±11.7 



lb: rr = 164.4 ±8.1 



Ila : Tc == 194.9 ± 7.4 



lib : Tc = 159.5 ± 7.8 



rl/3 



±5.2 



iMd 



±9.2 



iMd' 



,1/3 



±9.4 



(3.12) 
(3.13) 
(3.14) 
(3.15) 



where the uncertainties come from the quantities |g.i/3 and |^^ as denoted by the subindices, the possible 
variations of \„i/3 and |^^ arc given in TablctO and |31 

From the above analysis and results, we would like to address the following points: First, we may draw 
attention to the effect of the VEV v{z). From Fig[21 it is seen that the effects in all models are similar 
as long as it satisfies the boundary conditions (|2.6p and (|2.7p with the given a, regardless of the form of 
v{z). On the other hand, it turns out that the appearing of the peaks moves to high temperature as a^'^ 
increases. Next, let us turn to the effects of iii. From FiglH it is noticed that the behaviors of the peak 
become softer in all models, and the locations of the peak shift to the high temperature as /i^ increases. 



As for the effect of quark mass, considering the constituent quark mass defined via m^ 



a|^,■' 



witli iriq and cr / y? denoting the current quark mass and the dynamically generated quark mass via chiral 
condensation, respectively. As the relation ea. (j2.8|) is kept to be held, the current quark mass niq and 
the dynamically generated quark mass or chiral condensation a are varying oppositely, which can be seen 
explicitly from TablelSl It can also be seen from Fig |5] that when the constituent quark mass fhq is going 
to be increased; which is toward the increased a associated with the decreased mq due to the relation 



ea. (|2.8p . the behavior of the QNS is getting steeper, and vice versa. This is consistent with the general 
point of view that the more massive the current quark mass becomes, the more moderate the behavior 
of the peak for the QNS x/T^ gets. In this sense, it can be interpreted that the effect of current quark 
mass iTiq is more effective than the one of dynamically generated quark mass or the chiral condensation a. 
Finally, we plot in Fig|4]the results in all models with parameters given in Tableland the temperature is 
normalized by the critical temperature Tc, so that each peak denoting Tc come to 1 on the x-axis. 




Figure 4: The results of all models with the parameters given in Table[2l The curves corresponding to different 
models may be distinguished from their locations corresponding to large value of T/T^ > 1.5, form top to bottom, 
they are; lb (red, dashed), la (green, solid), lib (Magenta, dashed), Ila (blue, dot-dashed). The x-axis and y-axis 
denote the normahzed temperature T/Tc and x/T^ -, respectively. 



4 Summary and Remarks 

In this paper, we have extended the predictive soft-wall AdS/QCD model [13] to a finite temperature 
system and carried out the calculation for the quark number susceptibility with finite quark mass. It has 
been shown that the critical behaviors in all models agree well with the lattice QCD calculation [15,16]. 
This further demonstrates that the predictive soft-wall AdS/QCD model proposed in [13] is a reasonable 
and realistic effective model for describing the QCD confinement and chiral symmetry breaking, where 
the conformal symmetry breaking at infrared region due to nonperturbative QCD effects is characterized 
by the modified infrared behavior of metric parameterized by ^^ in the metric and the nontrivial dilaton 
background (namely the soft-wall) which is determined by both the IR parameter /i^ and the VEV v{z). 
It seems reasonable to suppose that v{z) has almost the same effects as long as it satisfies the boundary 



condition (j2.6p and (j2.7p . In addition, it has been noticed that the critical temperature is put to a low 
temperature by the soft-wall, while the effect of IR behavior /i^ will pull it up to high temperature. Thus 
the matching of the critical behaviors can be seen as the result of a very fine balance between /i^ and the 
soft- wall. As for the effect of quark mass given by the constituent quark mass rhq = iriq + ct//!^ , it has been 
shown that the effect of current quark mass nig is more significant than the one of chiral condensation a 
or dynamically generated quark mass when the relation cq. (|2.8p is imposed, which may be explicitly seen 
from the relation fhq = ruq + f^m'^/{2mq^'^). 

Here we only consider two fiavors, it would be interesting to investigate the nature of QCD transition 
with three flavor case, and also with including the chemical potential. On the other hand, the quartic 
interaction term of the scalar field has been dropped in eq. (|2.3p . while it has been shown that it has a 
sizable role in the computation of the mass spectra of resonance mesons [13]. Also the back reaction effects 
of the dilaton field and gravity may be the core in the improvement of our model [21]. We shall consider 
those interesting topics in our future studies. 
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